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ABSTRACT 

We provide a set of numerical N-body simulations for studying the formation of the outer Milky Ways's 
stellar halo through accretion events. After simulating minor mergers of prograde and retrograde 
orbiting satellite halo with a Dark Matter main halo, we analyze the signal left by satellite stars in 
the rotation velocity distribution. The aim is to explore the orbital conditions where a retrograde 
signal in the outer part of the halo can be obtained, in order to give a possible explanation of the 
observed rotational properties of the Milky Way stellar halo. Our results show that, for satellites more 
massive than r^ 1/40 of the main halo, the dynamical friction has a fundamental role in assembling 
the final velocity distributions resulting from different orbits and that retrograde satellites moving 
on low inclination orbits deposit more stars in the outer halo regions end therefore can produce the 
counter-rotating behavior observed in the outer Milky Way halo. 
Subject headings: galaxies: kinematic and dynamics; Galaxy: halo; methods: numerical 



1. INTRODUCTION 

The Galactic halo has long been considered a sin- 
gle component. However, evidences in the past few 
decades have i ndicated that it may be more complex 
(iPreston et al. . 1991: Maicwski 1992; Kinman ct alj 
19941: iCarneyetall Il996l: iChiba fc Beers I IMJol : 
Kinman et al. I 120071 : iMiceU et all I2008D . Recently 



Carollo et al. I ()2007[ ) confirmed the existence of a two- 



components halo, analyzing a large sample of calibration 
stars from the Sloa n Digital Sky Survey (SDSS) DR5 
(jAdelman-McCartv et al. II2007D . According to Carollo 
et al. the Galactic halo consists of two overlapping 
structural components, an inner and an outer halo. 
These components exhibits different spatial density pro- 
files, stellar orbits and stellar metalicities. In particular 
the inner-halo stars show a small (or zero) net prograde 
rotation around the center of the Galaxy. Outer halo 
stars possess a clear retrograde net rotation. 

The theory of galaxy formation in a Lambda cold 
dark matter (LCDM) Universe predicts galactic stellar 
halos to be built from multiple accreti on events start- 
ing from the first structure s to collapse (jWhite fc Rees I 
119781: iSearlefc Zi^ll978D . Halos, composed of both 
dark and baryonic matter, grow by merging with other 
halos. While the gas from mergers and accretions loses 
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its energy through cooling and settles into a disk, the 
non dissipative material (accreted stars and dark matter) 
form a halo around the Galaxy. In a LCDM universe, a 
dark matter halo big enoug h to host the Milky Wa y con- 
tains 300-h/-100 subhalos (jDiemand et al. 1120041) (such 
a large number of subhalos is rela ted to the s o-called 
"missing satellites problem" . iMoore et al. iri999[ ). 

There is significa nt evidence of past accretions into the 
Milky Way (.Yannv"eran l2003t iBelokurov et afl 120071: 
iGrillmair II2OO90 . and therefore we will accept the sce- 
nario that the Milky Way halo has been formed by 
subhalos accretion and we focus on possible origin of 
a retrograde (counter-rotating) outer stellar halo, de- 
fined as the set of stars orbiting around the Galaxy 
at large distance fro m the disk plane, say z > 15 kpc: 
I Carollo et al. I (|2007 D found that this distance separates 
the inner and the outer halo. It is expected that such 
halo is made of stars strippe d from me r ging o r disrupted 
satellites of low mass (Zolotov et al. I (|2009D : for a de- 
tailed study of accreting events onto a Mi lky- Way like 
halo, see also iBovlan-Kolchin et al. 1 ()2010[ )). However, 
from cosmological simulation, prograde an d retrograde 
mergers are approxima tely equally likely (jSales et al. I 
[2007tlRead eTaI~l l200g). 

In order to determine if, and in what condition, we 
can obtain a retrograde signal in the stellar distribution 
in the outer halo, we simulate minor mergers of a satellite 
halo onto a main halo of Galactic mass, with a mass ra- 
tio Mprimary/AfsatoUitc ~ 40. Wc put the Satellite on two 
orbits, one prograde and one retrograde. After the satel- 
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lite completes its merging/disruption process, we identify 
stars at distance larger than 15 kpc from the disk plane 
and simply analyze their rotation velocities. Our main 
finding is that a counter-rotating stellar halo naturally 
arises when minor mergers happen on orbits with a low 
inclination with respect to the disk plane. 

The plan of the Letter is the following. In Section[51 we 
describe our simulations; in Section[3l we give our results 
on counter- rotating outer halo stars, and in Section 2] we 
draw our conclusions. 

2. SIMULATIONS 

We use a primary Dark Matter (DM) halo containing 
a stell ar, rotating exponent ial disk. The DM halo has a 
NFW ([Navarro et al. 1 0-997') radial density profile, and a 
mass, radius and concentration appropriate for a Milky 
Way like DM halo at redshift z — 0. DM particles have 
velocities given b y the local equilibrium approximation 
(|Hernquist II1993D . Into our halo, we embed a truncated 
stellar disk, having an exponential surface density law: 
Pstars = Poexp— (r/ro) where ro is the disk scale length, 
and po is the surface central density. We obtain each 
disk particle ' s pos ition using the rejection method by 
iPress et al. I ()1986[ ): the dis k is in gravitational equ ilib- 
rium with the DM halo (see ICurir fc Mazzei I (|1999D for 
further details). We choose for the minor merger satel- 
lite a mass ratio of ~ 40, similar to the estimated mass 
ratio of the LMC to the Milky Way halo. The satellite 
contains a stellar bulge, with a Hernquist radial density 
profile. We realized our DM-|- bulgc satellite configura- 
tion as in IVillalobos fc Helmll (J2008). All the physical 
parameters of our merger are listed in Table [T] 

We simulated prograde mergers, in which a satellite co- 
rotates with respect to the disk spin, and retrograde ones 
with a count er-rotat i ng satellite. We chose two orbits 
used in Read et al. I (|2008l for studying the thickening 
of the disk due to the same kind of minor merger: a 
low-inclination one, with a 10 degree angle with the disk 
plane, and a high inclination one with a 60 degree angle. 
Initially, the center of the primary halo stays in the origin 
of our coordinate system and the satellite is in (x,y,z) = 
(80.0, 0.27, -15.2) kpc for the low-inclination orbit and 
(15.0, 0.12, -26.0) kpc for the high-inclination one. The 
(x,y,z) components of the velocity of the satellite are, 
in the prograde case, ( 6.3, -62.5, 0.35) km/s for the 
low-inclination orbit and (-1.2, 80.1, 2.0) km/s for the 
high-inclination one. The retrograde orbits have the y- 
component of the velocity inverted. 

The secondary has always a spin parameter A = 0, 
where we define A = J/{^2MVR), with M being the 
mass inside a radius R and V = GM/R the circular ve- 
locity, as in iBullock et al. I (J200ll ) . We use simulations 
with a primary halo with a spin parameter A = 1 (simu- 
lations "A" ) and A = (simulations "B" ) . We assign the 
angular momentum to DM particles using a rigid body 
rotation profile. The angular momentum of DM particles 
is always aligned with that of the stellar disk. 

Our primary halo has 10^ DM particles inside the virial 
radius (^ 2.5 • 10® in total) and one million star parti- 
cles in the exponential disk, with mass Mum = 10® M© 
and M, = 5.97 x lO'^ M© respectively. Our DM-|-Bulge 
secondary halo has 1.1 • 10^ DM particles and 10*^ bulge 
star particles, with masses Mgat = 1-95 x 10^ M© and 

Mhnl.n = 2.38 X 10^ Mr; 



^ bulge 



i©. 



We use a Plummer-equivalent gravitational softening 
length e — 0.5 kpc, and e = 0.25 kpc for bulge star 
particles. We run all our si mulation using the public 
parallel Treecode GADGET2 (iSoringel I I2005I). 

To test the convergence of our results with resolution, 
we re-run our set "A" with ten times more particles in 
the satellite halo. We also changed the spin parameter of 
primary and secondary h alos, its radial distrib ution using 
the profile from Bullock (jBullock et al. 11200 It ) , and their 
coupling. 

We analyze positions and kinematic of the bulge stars, 
once the satellite completes it merging with the primary 
halo ^. We run all our simulations for t — 4.63 Gyrs, 
corresponding to ~ 16 dynamical timescales of the main 
halo. 

3. RESULTS 

In all cases, in our simulations the satellite is slowed 
down by dynamical friction exerted on it by both disk 
and halo particles. At the first pericenter, tidal forces 
deform the satellite, redistributing its particle's energy 
(violent relaxation), and strips away some stellar parti- 
cles (tidal stripping). The process continues at each sub- 
sequent passage to the pericenter, until no recognizable 
self-gravitating structure is present anymore. Stripped 
star particles tracks the orbital pattern of the satellite. 
The spatial distribution of stars depends upon the dy- 
namical history of the satellite to which they belonged, 
i.e., how quickly it loses its orbital energy, how strongly 
it gets disrupted by tidal forces, and how these effects 
modify the orbit itself during the merger process. 

We rotate our coordinate system so that the stellar disk 
lies in the X-Y plane. The origin is in the disk center of 
mass. To detect a rotation signal in the outer halo stars, 
we simply took all the star particles initially belonging 
to our satellite, and having a coordinate Z > 15 kpc or 
Z < —15 kpc. We then calculate the rotation velocity of 
such particles in the disk plane. 

3.1. Distributions of Vtf, 

In Figure [l] we show histograms of the rotation ve- 
locity obtained in the four simulations of the set A: the 
satellite is either co-rotating or counter-rotating with re- 
spect to the disk, and the orbit has either a low or a 
high inclination with respect to the disk plane. In the 
upper panels, we show histograms at the final time of 
our simulations; in the lower ones, we performed an av- 
erage over five consecutive snapshots, (^ 100 Myr) to 
get rid of possible sampling effect on the particle orbits 
which can rise using a particular instant of time. We 
show both the distribution of stars from single simula- 
tions, and the one resulting from taking together star 
particles from prograde and retrograde orbits having the 
same initial inclination. The latter gives us a hint on the 
possible rotation signal origin in the case a similar num- 
ber of prograde and retrograde minor accretion events 
happen during the formation history of a galaxy, under 
the hypothesis that their inclination is the same. It is 
immediately clear from the Figure that our couple of low 
inclination orbits show an excess of counter-rotating stars 
in the outer halo. 

^ We define a merger to be complete when the z coordinate of 
the CM of all satellite stars remains within 2 kpc from the initial 
disk plane. 
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Table 1 
Properties of main and satellite halos. Column 1: Halo. Column 2: Virial mass, in M©; here, we refer all 

OUR virial quantities to an OVERDENSITY of 200 TIMES THE MEAN COSMIC DENSITY. COLUMN 3: ViRIAL RADIUS, IN 

KPC. Column 4: NFW concentration parameter. Column 5: Truncation radius of the main halo; the 

SECONDARY HALO HAS AN EXPONENTIAL CUT-OFF IN DENSITY. COLUMN 6: DiSK SCALE RADIUS (mAIN HALO ONLY) . COLUMN 

7: Disk truncation radius, in kpc. Column 8: Disk stellar mass, in Mq. Column 9: Hernquist scale radius, in 

KPC. Column 10: Mass of the stellar bulge, in Mq. 
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Fig. 1. — Histograms of rotation velocities for star particles in 
the outer halo, at the end of the set A of merger simulations. Upper 
panels show the histograms at the final time, lower panel show the 
same histograms, averaged over five simulation snapshots. In red 
(dashed line), we plot the histogram obtained for our retrograde 
orbits; in green (dotted line), those for our prograde orbits. In 
black (solid line) we show the sum of the two. Left column is for 
the low inclination orbits, right column for the high inclination 
ones. We used 50 velocity class, equispaced, between -300 and 300 
km/s in all histograms. 
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Fig. 2. — Upper row: histograms of rotation velocities for star 
particles in the outer halo, at the end of the set B of merger simu- 
lations, in which the main halo has no spin. In red (dashed line), 
we plot the histogram obtained for our retrograde orbits; in green 
(dotted line), those for our prograde orbits. In black (solid line) 
we show the sum of the two. Left panel is for the low inclination 
orbits, right panel for the high inclination ones. We also show the 
histogram for the case of spin 1 (dashed-dotted lines). Lower row: 
histogram of rotation velocities for star particles in the outer halo, 
for our set A, at the basic resolution ("LR", black continuous line) 
and at a ten times better mass resolution ("HR", red dashed line). 



Table [5] reports the number of satellite star particles in 
the positive and negative peak in all of our cases, and the 
number of particles having rotation velocity frot < — 10 
km/s and Vmt > 10 km/s. It also reports the same num- 
bers obtained from our set B, in which the DM halo has 
no spin. We give such numbers for the final time of our 
simulations. The counter-rotating excess signal, defined 
as the fraction of the number of counter- to co- rotat- 
ing satellite stellar particles, is 3.39 at peaks and 1.98 
overall for the low- inclination orbits, while it is 1.3 (at 
peaks) for the high-inclination ones. Note that, for high- 
inclination orbits, the retrograde case shows a peak at a 
rotation velocity near to Wrot = 0. This is because high 
inclination orbits have a larger impact parameters. The 
disk responds by tilting more than in the low-inclination 
case, and as a result, rotation velocities in the disk refer- 
ence frame are shifted towards more positive values. This 
effect is also present in our low-inclination runs, but it is 
much smaller because the tilt of the disk is smaller. 

It remains to be determined if the counter-rotating sig- 
nal is caused by the interaction with the stellar disk of 



the main halo, or by the DM of the halo itself. 

The upper row of Figure [5] shows the results of the 
same analysis performed on our set A, but in the case in 
which no halo spin is present (set B). Here we only show 
the results for our last snapshot; as in Figure [3 averag- 
ing over five snapshots makes no appreciable difference. 
Also in our set B, low inclination retrograde orbits pro- 
duce more counter-rotating star particle than co-rotating 
stars produced by prograde orbits. Again, this is not the 
case for high inclination orbits. From Figure [2l it is clear 
that the excess of counter-rotating star is clearly smaller: 
and this is due to the fact that more co-rotating stars are 
produced by prograde orbit in our no spin case than in 
spin 1 case. From Table [21 the counter-rotating excess 
signal is 1.60 at peaks and 1.59 overall in the low inclina- 
tion case, and 1.17 at peaks and 1.02 overall in the high 
inclination one. 

Therefore, both disk rotation and halo spin contribute 
to the slowing-down of prograde orbits and to the con- 
sequent smaller amount of high-energy star particles 
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Table 2 
Number of outer halo satellite star particles in the positive and negative peak in all of our cases, and the 

NUMBER OF OUTER HALO SATELLITE STAR PARTICLES HAVING ROTATION VELOCITY Vrot < — 10 KM/S AND lirot > 10 KM/S IN 

OUR LOW INCLINATION RUN. RESULTS ARE FOR THE FINAL TIME OF OUR SIMULATIONS. FOR EACH INCLINATION, WE 

CONSIDERED TOGETHER PROGRADE AND RETROGRADE MERGERS. FiRST ROW: NUMBERS OF SATELLITE STAR PARTICLES 

HAVING ROTATION VELOCITIES SMALLER THAN —10 KM/S OR HIGHER THAN 10 KM/S, FOR THE LOW INCLINATION CASE. 

Second row: number of star particle in the positive and negative peak, for the low inclination case. Third 
row: number of star particles in the positive and negative peak, for the high inclination case. second 

column: number of counter- ROTATING SATELLITE STARS WHEN THE DM HALO HAS A SPIN A = 1. THIRD COLUMN: 

NUMBER OF CO-ROTATING SATELLITE STARS WHEN THE DM HALO HAS A SPIN A = 1. FOURTH COLUMN: NUMBER OF 

COUNTER-ROTATING SATELLITE STARS WHEN THE DM HALO HAS NO SPIN, A = 0. FiFTH COLUMN: NUMBER OF CO-ROTATING 

SATELLITE STARS IN THE SAME A = CASES. 



Distribution 



A = 1 counter A = 1 co A = counter A = co 



lOdeg TOTAL v, > |10| km/s 20292 10237 

lOdeg PEAK 6510 (-54) 1917 (-1-54) 
60deg PEAK 3581 (-1-6) 4010 (64) 



16670 13033 

4455 (-54) 3419 (54) 
2741 (-6) 3071 (78) 



stripped from satellites that can reach the outer halo. 

The lower row of Figure [2] shows that our result does 
not depend on the mass resolution of our satellite halo. 
Even using 10 times more particles in the secondary, only 
our low inclination couple of mergers shows an excess of 
counter-rotating stars in the outer halo. 

3.2. Description of the mergers 

Two effects are acting on the satellite halo in these 
kind of mergers: dynamical friction and tidal disrup- 
tion. The first one is exerted both by the main halo 
DM particles and by the disk star particles. The sec- 
ond is most important near the center of the main halo, 
where the gravitational potential is stronger. But, dy- 
namical friction depends upon the details of the satel- 
litc o r bits. It is already k nown (Quinn fc Goodma _n1 
lOM iWalker et al. I fl996t iHuang fc CarlbergI Il997[ 
Velazquez fc White I [1991 IViUalobos fc Helmi I 120081) 



that prograde orbits tends to decay faster than retro- 
grade ones. The dynamical friction force goes as i^dyn oc 
1/v^ (iBinnev fc Tremaine 112008 ) . where Vs is the velocity 
of the satellite relative to the field particles; retrograde 
satellites have higher Vs with respect to prograde one, 
since in the first case the rotation velocity of the satellite 
is opposite to that of the disk and the main DM halo 
particles. Particles stripped from the satellites will re- 
main on the orbit on which the satellite was when they 
were stripped. Therefore, retrograde satellites "deposit" 
more star particles in the outer halo regions, producing 
the signal we observe in our simulations, since their or- 
bits have a longer decay time. Obviously, to obtain this 
effect, the orbital velocity must lie on the disk plane ^, 
which is almost the case for our low inclination orbits. 
High inclination encounters do not show such a behavior. 

In Figure [3l we plot the position of the center of mass 
of all the star particles belonging to the satellite as a 
function of time. In our low inclination mergers (left 
panel), the prograde orbit decades faster: already after 
1.5 Gyr, the difference with the behavior of the retro- 
grade satellite is clear. In the high inclination cases the 
effect, though present, is much smaller, as expected. 

We also verified that the effect we show here does not 
depend upon the exact angular momentum distribution 

^ DM spin is parallel to the disk spin: also the effect of the DM 
halo spin is maximum in such a case 
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Fig. 3. — Position of the center of mass of all star particles 
belonging to the satellite, as a function of time. Left panel shows 
the result for our low inclination orbits, right panel for the high 
inclination ones. The center of mass position makes sense only 
until the satellite is not disrupted; this happens after ~ 3 Gyr 
in our low-inclination prograde simulation, ^ 4 Gyr in our low- 
inclination retrograde one, while in the high inclination cases cases 
the satellite is disrupted already after ~ 1 Gyr. 



of particles in DM halo, by repeating our experiments 
using a Bullock (2001) profile instead of a rigid rotation 
one for the DM rotation velocity. Finally, the initial dis- 
tance at which we put our satellite does not have a great 
impact on our results: we repeated the low-inclination 
mergers, putting the satellite at coordinates (29.5, 0.27, 
-5.2) kpc, thus much nearer to the center, and obtained 
the same qualitative behavior. 

4. CONCLUSIONS 

We performed controlled numerical simulations of mi- 
nor merger events. Our aim was to determine if an excess 
of counter-rotating stars in the outer stellar halo can be 
produced by couples of mergers having orbits with identi- 
cal inclination and opposite initial rotation. Our primary 
halo has an exponential stellar disk and its dark matter 
has a NFW radial density profile; our satellite's DM has 
an NFW density profile, and contains a stellar spheroid 
with a Hernquist density profile. The mass ratio of our 
merger is Mprimary/Msatciiitc ~ 40. We simulate low in- 
clination encounters, in which the angle of the satellite's 
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orbit with the disk plane is 10 degrees, and high inclina- 
tion ones (60 degrees). We use prograde and retrograde 
orbits and also vary the DM spin parameter of the main 
halo. DM spin is aligned with the disk rotation. Here 
we define all stars having a distance 2 > 15 kpc from the 
disk plane as belonging to the outer stellar halo. 
Our main results are the following: 

• low inclination mergers do produce an excess of 
counter-rotating satellite stellar particles in the 
outer halo, independently on the spin of the DM; 

• our 60 degrees mergers with our 1:40 mass ratio 
are not able to give a significant counter-rotating 
signal, also owing to the disk tilt produced by the 
merger itself; 

• the fraction of counter- rotating to co- rotating satel- 
lite stars in the outer stellar halo is higher if the DM 
has spin. 

From our controlled experiments, we now have an in- 
dication of the possible origin of the counter-rotation of 
stars observed in the Galactic outer halo in the context 
of a hierarchical ACDM model of galaxy formation, in 
which many major mergers happen during the history 
of the Galaxy. Even if, statistically, the number of pro- 
grade and retrograde minor merger is the same, still a 
counter-rotating signal can arise if such mergers predom- 
inantly happen along low inclination orbits. Since mat- 
ter accretion, in a CDM dominated Universe, mainly oc- 
curs along filaments, this will be the case if the galaxy 



disk is co-planar to the (majority of) filaments. The 
disk-filament alignni ent issue is still debated (see e.g. 
iBrunino et al. I [20071 ): from our results, we expect that 
if the galactic disk were perpendicular to the main ac- 
cretion streams, no counter-rotating signal should be ob- 
served in the outer halo star distribution. 

Our orbits are not cosmologically motivated, since the 
aim of our experiment is to determine if and in what cases 
an excess of counter-rotating stars in the outer halo can 
be produced. Of course, in a realistic case mergers will 
occur with orbits having a number of different inclina- 
tions, giving rise to a velocity distribution which will not 
show such a clear, double peaked signal as that detected 
in the present work. A detailed study of the orbital pa- 
rameters of minor mergers in a statistically significant 
set of cosmological Galaxy-sized halo cosmological accre- 
tion histories is needed before a quantitative comparison 
between theory expectations and observation can be per- 
formed. 



We acknowledge useful discussions with Daniela Car- 
oUo, Stefano Borgani, Gabriella De Lucia, Antonaldo Di- 
aferio and Alessandro Spagna The simulations were car- 
ried out at CASPUR, with CPU time assigned with the 
"Standard HPC grant 2009" call, and at the "Centro In- 
teruniversitario del Nord-Est per il Calcolo Elettronico" 
(CINECA, Bologna), with CPU time assigned under an 
INAF/CINECA grant. A.C. acknowledge the financial 
support of INAF through the PRIN 2007 grant n. CRA 
1.06.10.04 "The local route to galaxy formation". 



REFERENCES 



Adelman-McCarthy, J. K. et al, 2007, ApJS, 172, 634 

Belokurov, V. et al, 2006, ApJ, 642, L137 

Binney, J., & Tremaine, S, 2008, Galactic Dynamics: Second 

Edition, Princeton University Press, Princeton 
Boylan-Kolchin, M., Springel, V., White, S. D. M., Jenkins, A., 

2009, arXiv/091 1.4484 
Brunino, R., Trujillo, I., Pearce, F. R., & Thomas, P. A., 2007, 

MNRAS, 375, 184 
Bullock, J. S., Dekel, A., Kolatt, T. S., Kravtsov, A. V., Klypin, 

A. A., Porciani, C, &Primack, J. R., 2001, ApJ, 555, 240 
Carney, B. W., Laird, J. B., Latham, D. W., &AguiIar, L. A., 1996, 

AJ, 112, 668 
Carollo, D. et al., 2007, Nature, 450, 1020 
Chiba, M., & Beers, T. C, 2000, AJ, 119, 2843 
Curir, A., & Mazzei, P., 1999, A&A, 352, 103 
Diemand, J., Moore, B., &StadeI, J., 2004, MNRAS, 352, 535 
Grillmair, C. J., 2009, ApJ, 693, 1118 
Hernquist, L., 1993, ApJS, 86, 389 
Huang, S., feCarlberg, R. C, 1997, ApJ, 480, 503 
Kinman, T. D., Suntzeff, N. B., &;Kraft, R. P., 1994, AJ, 108, 1722 
Kinman, T. D., Cacciari, C, Bragaglia, A., Buzzoni, A., ifeSpagna, 

A. 
Majewski, S. R., 1992, ApJS, 78, 87 
Miceh, A. et al., 2008, ApJ, 678, 865 
Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, 

J., & Tozzi, P, 1999, ApJ, 524, L19 



Navarro, J. F., Frenk, C. S., & White, S. D. M. 

Press, W. H., Flannery, B. P., &Teukolsky, S. A., 1986, Numerical 

recipes. The art of scientific computing, Cambridge: University 

Press, 1986 
Preston, G. W., Shectman, S. A., & Beers, T. C, 1991, ApJ, 375, 

121 
Quinn, P. J., & Goodman, J., 1986, ApJ, 309, 472 
Read, J. I., Lake, G., Agertz, O., & Debattista, V. P., 2008, 

MNRAS, 389, 1041 
Sales, L. V., Navarro, J. F.; Abadi, M. G., Steinmetz, M., 2007, 

MNRAS, 379, 1464 
Searle, L., & Zinn, R., 1978, ApJ, 225, 357 
Springel, V., 2005, MNRAS, 364, 1105 

Velazquez, H., &; White, S. D. M., 1999, MNRAS, 304, 254 
Villalobos, A., & Helmi, A., 2008, MNRAS, 391, 1806 
Walker, I. R., Mihos, J. C, & Hernquist, L., 1996, ApJ, 460, 121 
White, S. D. M., & Rees, M. J., 1978, MNRAS, 183, 341 
Yanny, B. et al., 2003, ApJ, 588, 824 
Zolotov, A., Willman, B., Brooks, A. M., Governato, F., Brook, 

C. B., Hogg, D. W., Quinn, T., & Stinson, G., 2009, ApJ, 702, 

1058 



